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Tunable spin-polaron state in a singly clamped semiconducting carbon nanotube 
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We consider a semiconducting carbon nanotube (CNT) laying on a ferromagnetic insulating sub¬ 
strate with one end depassing the substrate and suspended over a metallic gate. We assume that 
the polarised substrate induces an exchange interaction acting as a local magnetic field for the elec¬ 
trons in the non-suspended CNT side. Generalizing the approach of I. Snyman and Yu.V. Nazarov 
[Phys. Rev. Lett. 108 , 076805 (2012)] we show that one can generate electrostatically a tun¬ 
able spin-polarized polaronic state localized at the bending end of the CNT. We argue that at low 
temperatures manipulation and detection of the localised quantum spin state is possible. 

PACS numbers: 85.85.+j,71.38.-k, 85.75.-d 


Nanoelectromechanics with suspended carbon nan¬ 
otubes evolved rapidly in last few years m • Recently 
I. Snyman and Yu.V. Nazarov [8] considered a semicon¬ 
ducting CNT laying on an insulating substrate with one 
end of it suspended. A metallic gate below both the in¬ 
sulating substrate and the suspended part of the CNT 
generates an homogeneous electric field (cf. Fig. 1 of Ref. 
0 )- The mechanical bending of the suspended part of 
the nanotube induces then a spatial inhomogeneity of 
the electrostatic potential along the CNT forming a min¬ 
imum at the deformable end of the wire. The competition 
between such an electrostatic bending with both the elas¬ 
tic potential of the CNT and the quantum rigidity of the 
electronic wave function makes the mechanical bending 
as well as the formation of the localized polaronic state 
at the movable end of the CNT to occur as a first order 
phase transition as a function of the electric field. The 
estimate for the critical field for a realistic experimental 
set up was predicted in Ref. [8] to be 0.01 V//im. 

An impressive effort of the nanoelectronics community 
is currently deployed to manipulate and exploit the elec¬ 
tronic spin degrees of freedom in transport devices (spin- 
tronics) [9]. In this context the possibility of magnetic 
gating , i.e. the use of ferromagnetic leads inducing mag¬ 
netic exchange fields F ex / fi b (with /i# the Bohr magne¬ 
ton) on the electronic spin is currently actively investi¬ 
gated [TUHS]. More surprisingly such exchange fields can 
have remarkable consequences also on the dynamics of a 
nano-mechanical system for which dynamical (shuttle) 
instabilities, strong spin-polarized currents, and cooling 
have been predicted [T411T6] 

In this paper we show that the system discussed by 
Snyman and Nazarov |8] in presence of a magnetic di¬ 
electric substrate allows the formation of a localized fully 
polarized polaronic state. The experimentally observed 
exchange energy E ex (see Ref. [131 14.) turns out to be 
as large as tens of Kelvins, thus being of the same order 
of magnitude of the localization energy for an electron in 
a CNT on the scale of the micrometer. This allows for 
a high tunability of the polaronic state by means of two 
electric gates, below the suspended and non-suspended 
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FIG. 1. Schematic of the system considered: a CNT laying 
on a magnetic substrate ( C ) and protruding out of a length 
L. Two independently adjustable gates (Vbi and Vg 2 ) are 
shown and a contact ( C ) on the substrate side. The potential 
for spin up and down ([/+ and U- ) is also sketched above. 


part of the CNT (see Fig. [I]). As a result a continous 
electrostatic tuning of the localization length and the 
bound state energy can be achieved, forming a stabil¬ 
ity diagram of spin-up and spin-down polaronic states. 
Detection of the state of the system can be envisaged by 
use of a nearby single-electron transistor, for which the 
CNT tip acts as a gate m- Fully electric manipulation 
of the mechanical and electronic spin state of the CNT 
is thus possible in this system. 

The system. Following Ref. [8j let us consider a CNT 
laying on a substrate with a suspended part protruding 
out of a length L (see Fig. 1). In Ref. [8] it has been 
shown that the wavefunction ^(x) of the electrons in the 
valence band of the CNT can be describled by a standard 
one-dimensional Schroedinger equation with an effective 
mass m* = 0.6ra e ao/r, where m e is the electronic mass, 
ao the Bohr radius, and r the radius of the CNT. The 
variable x parametrizes the position along the CNT, its 
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value is 0 at the edge of the substrate and L on the tip 
of the suspended part. The length of the CNT on the 
substrate is supposed to be L, and is taken infinity 
for simplicity. Then vanishing boundary conditions ap¬ 
ply at x = L and x = — oo. As in Ref. [8] the CNT 
can bend with a displacement y{x) (for 0 < x < L) in 
the direction orthogonal to the substrate. The elastic en¬ 
ergy cost reads IY f dx [y"(x)] 2 / 2, where I = 6.47raor 2 
is the second moment of area of the tube cross section 
and Y the CNT Young modulus (of the order of 1.2 
TPa). Single clamping implies that y(Q) = y'( 0) = 0 
and y"(L) = y"'(L ) = 0. In this paper we will restrain 
to the classical description of the deflection. The tip of 
the CNT on the substrate is in tunneling contact with a 
metal whose chemical potential can be tuned close to the 
valence band of the CNT by adjusting the electric poten¬ 
tial. Up to now the description followed closely Ref. [8]. 
We introduce now the main difference: We will assume 
that the substrate is a magnetic insulator that induces an 
exchange interaction term -E ex dxa\'ip a (x) | 2 for the 
electrons being in the CNT over the substrate {x < 0). 
The variable a indicates the spin projection in the z di¬ 
rection. This creates a spin-dependent potential so that 
the spin-up electrons (a = +) are attracted in the x < 0 
region. In order to tune the potential we assume that 
two different gates are present, one below the magnetic 
substrate and an other one under the suspended part. By 
changing independently the potentials on the two gates it 
is possible to modify the electrostatic potential V and the 
electric field E acting on the electrons on the suspended 
part (taking the non-suspended region as a reference for 
the potential, cf. Fig. [l]). We can then write the full 
Hamiltonian for the problem as follows (Q x is the Heavi¬ 
side function): 
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then diverges at the threshold for its appearance. We will 
thus see that the presence of V and F ex allows to create 
a spin-dependent tunable bound state, that is associated 
to a displacement of the CNT tip. 

Electronic problem. Let us begin with the purely elec¬ 
tronic problem [y(x) = 0 for all x\. The ground state can 
be found by solving the Schroedinger equation: 
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for each spin projection. The presence of a bound state 
is signaled by the existance of a solution of Eq. © with 
e a < —crE’ex the bottom of the relative band. Taking 
6(j < — crE’ex as a reference in energy the problem for 
each spin species reduces to that describled by Eq. © 
with E e x —> 0 and eV —eV — crE eyi = U. The solution 
can then be found by matching the wave function ^{x) = 
Ae KX for x < 0 with ip(x) = Be lkx -f Ce~ lkx for x > 0 
at x = 0 asking the continuity of the wave function and 
of its derivative. The boundary conditions lead to the 
eigenvalue equation e~ 2lkL = —{ik + E)/(ik — n) with 
kL = [—2m*e^/^ 2 ] 1 / 2 , kL = [2m* (U + e^/h 2 ] 1 / 2 and 
Cb < 0 the bound state energy. At the threshold —>• 
0 - , thus there k vanishes and the eigenvalue equation 
reduces to e~ 2zkL = — 1. This gives kL = 7r/2 and the 
threshold value for the potential U t = (■ tt/2) 2 Ek , with 
Ek = h 2 /(2m*L 2 ) the kinetic energy scale. For U — 
U t E k the localization length ft -1 = 2 LE K /(U — U t ) 
and diverges as anticipated. By changing U it is then 
possible to adjust the spread of the wave function on 
the magnetic substrate. Since the two spin species feel 
a different potential only on the substrate, this allows to 
change continuously the energy difference of the up and 
down bound states. The bound state energy for each spin 
state reads: 


e<r = -vE ex + e b (eV - aE ex ) (3) 


- ( E ex cr6- x + eV9 x - eEy(x )) |Vv0*0| 2 


( 1 ) 


The first term in Eq. 0 gives the quantum kinetic en¬ 
ergy, the second the elastic energy, and the third is a sum 
of three parts: the exchange energy, the electrostatic po¬ 
tential and its variation induced by the deflection [y{x)\ 
of the CNT. In Ref. |5], for V = 0 and for E ex = 0, it has 
been shown that it exists a critical value of the electric 
field E c for which the ground state is an electronic local¬ 
ized state on the CNT suspended part. The formation 
of the bound state is a first order transition: the CNT 
starts to bend only for E > E c and a metastable bound 
state exists for E c \ < E < E c . At E = E c the local¬ 
ization length is thus finite and typically much shorter 
than L. In order to have a tunable bound state it is nec¬ 
essary to have a smooth transition from the delocalized 
to the localized state. This is actually the typical case 
in quantum mechanics, by decreasing the depth of a po¬ 
tential well that allows a bound state one can delocalize 
progressively the wave function. The bound-state radius 


with the threshold value for V: eV a = U t + crE ex . A typi¬ 
cal picture of the eV dependence of the two bound states 
for E = 0 is shown dashed in Fig. [2] For U_ < V < U+ 
a unique bound state exists for spin down. Let’s de¬ 
fine V c the value for which the down-spin energy crosses 
the up-spin bottom band: e_(U = V c ) = —E ex . For 
U + < V < V c two bound states exist, but only the lowest 
one (spin up) is stable, since the spin-down lays above the 
bottom of the spin-up band, and any spin-flip perturba¬ 
tion allows its decay into the spin-up continuum. Finally 
for V > V c two stable bound states exist. Their energy 
splitting has a maximum at V c and then monotonically 
decreases as a function of V. This is due to the reduc¬ 
tion of the localization length reducing the effect of the 
exchange interaction that acts only for x < 0. Although 
both spin-up and spin-down polaronic states are stable 
at V > V c only one of them can be occupied due to the 
Coulomb blockade, whose repulsion energy turns out to 
be much larger than the polaronic bound state energy 
(^ Ek) at L 1 nm. This fact allows the formation of 
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FIG. 2. Dependence on the gate voltage V of the two spin- 
projection bound-state energy (e a ) for E^/Ek = 0.5, a = 0 
(i.e. E = 0, dashed) and a = 50 (full line). The bottom of 
the lower electronic band (spin-up) is shown dotted. 
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FIG. 3. For E ex = Ek and a = 175 gate voltage dependence 
of the bound state energies (top-left), fraction of localized 
charge or spin (top-right), deflection of the tip (bottom-left), 
and ratio of deflection to localized charge (bottom-right). 


a controllable single-electron fully spin polarized state at 
the protruding part of the CNT. 

Nanomechanical effects. We now consider how the sys¬ 
tem behaves when we let the CNT bend. It is not possi¬ 
ble any more to find the ground state energy analytically, 
we will thus follow closely the variational method used 
in Ref. |8| to which we refer for more details. We intro¬ 
duce the dimensionless variables z = x/L, h = H/Ek , 
/ = yYI/eEL 3 , <f> a = and the coupling param¬ 

eter a = {eE) 2 L 3 / (YIE x) • The problem can then be 
completely determined by giving only three independent 
coupling parameters: a, fi = E e ^/E Kl and v = eV/Ex • 
The functional to be minimized reads: 


h = 




(imt6- z + v8 z ) cpl 

(4) 


By writing <p{z) = Y^=i a niX ~ z ) n for z > 0 and 

J2n=i a n z n e KZ for z < 0, and f(z ) = J2n=i Kz n+1 one 
can enforce the boundary conditions and minimize nu¬ 
merically the functional in order to find the parameters 
{a n , b n , k} and thus the ground state energy e a with ex¬ 
plicit expressions for the bending and the wavefunction. 
The charge accumulated on the suspended part of the 
CNT in presence of an electric field induces a force that 
bends the tip of the CNT. The effective electronic poten¬ 
tial deepens and bending lowers the bound state energy. 
In particular it favors a stronger localization of the wave 
function on the tip (measured by ft -1 ). For E ex = eV = 0 
in Ref. [8] it is shown that the bound state forms with a 
first order phase transition for a > a c = 312.03. In or¬ 
der to keep a smooth transition we will consider the case 
a < a c and investigate the dependence of the bound state 
energy and wavefunction on eV/Ex for given values of 
E ex /Ex- 

Before considering the results of the numerical calcula¬ 
tion it is useful to estimate analytically the typical range 
of the displacement of the CNT tip induced by the lo¬ 
calization of the charge. Let us assume that the fraction 


n < 1 of an electron charge is accumulated on the CNT 
tip uniformely. A simple ansatz for the displacement is 
f{z ) = az 2 . It satisfies both the boundary conditions and 
the Euler equation f"" = 0. Substituing it in Eq. @ one 
has for the part proportional to a: 

h a = <a(2a 2 + n^~) (5) 

O 

this functional has a minimum at a = —n/12 = /(1). It 
gives a rough estimate of the dimensionless displacement 
of the tip by taking into account only the competition 
between the electric field and the elastic stiffness. The 
effect of the other two parameters is hidden into the value 
of n, that cannot be larger than 1. 

We present on Fig. [3] the numerical results for a = 175 
and fi = 1. One can see that the energy splitting of the 
two spin states is of the order of Ex = E ex (top-left 
panel). Defining n a = dzcf) 2 the fraction of charge 
(and spin) localized on the suspended part of the CNT 
one finds that for V = V c both bound states present a 
finite value of n a and n_ — n + « 0.17. The difference is 
slowly reduced for larger values of the gate voltage. The 
same can be said for the deflection of the tip of the CNT 
{fa = /(1) for each spin state, bottom-left panel). Fi¬ 
nally the bottom-right panel shows that the ratio f±/n± 
is actually close to the rough estimate 1/12. 

The plots of Fig. [3] show that a particularly important 
quantity is the value of the physical parameters (e CT , n a 
and f a ) at the threshold V c . The dependence on V is 
always monotonic and the maximum or minimum values 
are observed at V c . In view of manipulating the spin 
state, the value at V c gives thus a very good indication of 
the range in which the state can be accessible. We thus 
show in Fig. [4] as a function of a and for different values 
of (i the threshold V c , the energy splitting e_ — e+, the 
difference in the occupation n_ — n+, and in the deflection 
/_ — /+. As expected the critical voltage V c decreases as 
a function of <a, and in particular, for sufficiently small /q 
it vanishes when a approaches the critical value a c . The 
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FIG. 4. The four panels show as a function of a for /x = 
0.1,1., 10, 20 the critical gate voltage value V c (top left), the 
energy splitting of the two bound states e_ — e+ (top right), 
the difference in the fraction of localized charge n + —n- (bot¬ 
tom left), and the difference of the CNT tip deflection /_ —/+ 
(bottom-right); the last three quantities calculated at V = V c . 


bound-state energy splitting is monotone in a , since the 
electric field increases the localization of the bound state, 
and thus reduces the difference of the two states. Its a- 
dependence is rather weak. Even for /x ^ 1 the energy 
splitting remains of the order of that thus sets the 
main energy scale of the problem. Quite surprisingly the 
difference in the fraction of localized charge (n_ — n+) 
is not monotonic for small /x as a function of the electric 
field. This is due to the fact that the transition region is 
approached at different values of a for each spin state. A 
similar behavior is observed in /_—/+. One can conclude 
that the optimal value of a to observe a well defined 
bound state is between 100 and 200. 


Estimates. In order now to consider the possibility to 
observe the two bound states we discuss the typical scales 
of the problem. Expressing the radius in nm and the 
length in /xm Ek ~ 13.9(r/L 2 3 4 5 6 7 8 9 ) mK. The typical value of 
L ranges between 0.1 and 1/xm, leading to a range for Ek 
between few K to tens of mK, thus always accessible with 
standard cryogenics. The thermal and quantum fluctua¬ 
tions of the displacement of the tip plays also an impor¬ 
tant role, since they define the distinguibility of the dis¬ 
placement of the two bound states. From Eq. one can 
write an approximate potential for the tip displacement 
Sf — /(l) — /o (with /o = n/12 the equilibrium value): 
h a = 2 a(Sf) 2 . The equipartition theorem then gives 
for the thermal fluctuations Sfr = [F^T/(4<a£ 1 x)] 1 / 2 . 
Quantum fluctuations -S/q- has the same expression with 
ksT —* huOm. Since huim/Ex = 0.0332 independently of 
L or r ia mi then S/q = 0.09/v^- Expressing as above 
T in mK, L in /xm, and r in nm 5fr = 0.13I/[T/(r<a)] 1 / 2 . 
Those values have to be compared with /_ — /+ that 
are at best 0.04. /q is thus 5 times smaller of this value 
already for a = 100, while in order to keep 5fr small 
one needs T <C .09ra/L 2 . This is realizable for instance 
choosing L = .5/xm, r = 2nm, a = 200 and working at 
temperatures T ^ 20 mK (. Ek is 111 mK in this case). 

Conclusions. We have shown that combining elec¬ 
trostatic and magnetic gating the formation of a spin- 
polaronic state in a singly clamped CNT becomes possi¬ 
ble. Electric, magnetic, and mechanical tuning provides 
an effective manipulation of such spin-polaron states of¬ 
fering a controllable magneto-electro-mechanical trans¬ 
duction with single electronic charge and spin sensitivity 
involving sub nanometer mechanical displacement. 
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